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Abstract 

In this work we investigate the so-called model-space Brueckner-Hartree-Fock (MBHF) 
approach for nuclear matter as well as for neutron matter and the extension of this 
which includes the particle-particle and hole-hole (PPHH) diagrams. A central in- 
gredient in the model-space approach for nuclear matter is the boundary momentum 
beyond which the single-particle potential energy is set equal to zero. This is also 
the boundary of the model space within which the PPHH diagrams are calculated. 
It has been rather uncertain which value should be used for %m- We have carried 
out model-space nuclear matter and neutron matter calculations with and without 
PPHH diagrams for various choices of ku and using several modern nucleon-nucleon 
potentials. Our results exhibit a saturation region where the nuclear and neutron 
matter matter energies are quite stable as &m varies. The location of this region may 
serve to determine an "optimum" choice for kM- However, we find that the strength 
of the tensor force has a significant influence on binding energy variation with kM- 
The implications for nuclear and neutron matter calculations are discussed. 
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1 Introduction 



A common problem to non-relativistic nuclear matter calculations has been the 
simultaneous reproduction of both the binding energy per nucleoli (BE /A = 
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16 ± 1 MeV) and the saturation density p = 0.17 fm~ 3 . All modern realistic 
calculations of nuclear matter which employ two-body interactions such as 
Brueckner-like approaches [1-6], hyper-netted-chain [7] or coupled cluster [8] 
approaches are not able to reproduce in a satisfactory way the above nuclear 
matter data. If for a given interaction one is able to predict the saturation 
density, then the binding energy is underestimated. Similarly, if one is able to 
reproduce the binding energy of nuclear matter, then the saturation density 
is too high. This results in the so-called Coester band (see e.g. Fig. 3 below), 
where the saturation point for a given nucleon-nucleon (NN) interaction is 
correlated with the strength of the tensor force in the NN interaction. 

In the past years relativistic approaches to nuclear matter have also been de- 
veloped. For example, Brockmann and Machleidt [9,10] have recently shown, 
starting from a realistic NN interaction, that by performing a relativistic Dirac- 
Brueckner-Hartree-Fock (DBHF) nuclear matter calculation, one is able to 
meet the empirical nuclear matter data. Another line of approach is repre- 
sented by the relativistic mean field scheme of Serot and Walecka [11]. How- 
ever, the present relativistic approaches neglect the coupling to the negative 
energy solutions, and moreover relativistic approaches lead to a very small 
depletion of the Fermi sea [12]. This depletion is not consistent with recent 
(e, e'p) experiments. 

Therefore, we believe that it is still worth to properly investigate non-relativistic 
schemes based on extensions of the Brueckner-Hartree-Fock (BHF) approach. 
One of the main problems with Brueckner theory is that particle states are 
treated differently from hole states. In the standard BHF approach one em- 
ploys a discontinuous single particle (sp) spectrum, with a self-consistent BHF 
part for the nucleons below and a free-particle part above the Fermi surface. 
This leads to an unphysically large discontinuity of about 60 MeV for the sp 
spectrum at the Fermi level. 

Several methods have been proposed to overcome this difficulty Mahaux and 
his collaborators [6] used the real part of the reaction matrix to obtain a self- 
consistent sp spectrum which is continuous for all momenta. In the model- 
space Brueckner-Hartree-Fock (MBHF) approach proposed by Ma and Kuo[4] 
a continuous sp spectrum is obtained within a chosen model space (defined 
for momenta k < /cm), and the discontinuity at the model space boundary ku 
is rather small for proper values of ku- Although both methods yield more 
binding energy than the standard BHF approach, the nuclear matter empirical 
data are not reproduced. However, by including PPHH diagrams in the model 
space approach Song et al. [5] and Jiang et al. [13] have demonstrated that 
the saturation properties can be improved. Comparing with the standard hole- 
line expansion in BHF calculations, it is remarkable that the PPHH diagrams 
calculation can simultaneously increase the nuclear-matter binding energy and 
lower the saturation density p . 
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A central question concerning the above PPHH-diagram nuclear matter ap- 
proach is how to choose the boundary momentum ku- In the past [5,13], one 
usually took a "reasonable" value for the size of the model space such as ku 
= 3.2 fm _1 and proceeded to evaluate the nuclear matter PPHH diagrams. To 
our knowledge, PPHH-diagram nuclear matter calculations with other choices 
of ku have not been investigated in a systematic way. Are there certain criteria 
which may help determine the boundary k^p- Is there a saturation behavior 
in the sense that there exists a region where the results of PPHH-diagram 
nuclear matter calculations are insensitive to the choice of k M t The aim of 
this work is to study these questions. The answers to these questions may help 
confirm, or disconfirm, the rather encouraging results of the PPHH-diagram 
calculations reported earlier. 

The reader should note that another way of summing PPHH diagrams has 
been developed by Dickhoff, Polls and Ramos, see e.g., Refs. [1-3]. There a 
self-consistent Green's function (SCGF) approach is used, where one solves 
the Dyson equation in order to get the sp Green's function. Although the sp 
energies are kept real in the self-consistency scheme, as in our case as well, the 
two methods differ in the treatment of the poles in the in-medium scattering 
matrix. In the model-space approach, see the discussion in the two next sec- 
tions and Ref. [5], there are, by construction no poles at negative and positive 
starting energies in the effective interaction defined by the model-space reac- 
tion matrix. The poles from the retarded Green's function are circumvented by 
a contour integration [5]. One avoids thereby a numerically tedious evaluation 
of the two-body Green's function and effective interaction, as done in Refs. 
[1-3]. However, the SCGF method of Dickhoff, Polls and Ramos [1-3] allows 
one in a direct way to study other properties than just the binding energy. 
The one-body Green's functions can be used e.g., in the study of momentum 
distributions and spectral functions. 

Finally, our scheme to obtain an effective interaction appropriate for nuclear 
matter starts with a free nucleon-nucleon interaction V which is appropriate 
for nuclear physics at low and intermediate energies. At present there are sev- 
eral potentials available. The most recent versions of Machleidt and co-workers 
[14], the Nijmegen group [15] and the Argonne group [16] have a x 2 fit P er 
datum close to 1. The potential model of Ref. [14] is an extension of the one- 
boson-exchange models of the Bonn group [9], where mesons like 7r, p, i], 5, 
uj and the fictitious a meson are included. In the charge-dependent version of 
Ref. [14], the first five mesons have the same set of parameters for all partial 
waves, whereas the parameters of the a meson are allowed to vary. The re- 
cent Argonne potential [16] is a charge- dependent version of the Argonne V14 
[17] potential. The Argonne potential models are local potentials in coordi- 
nate space and include a 7r-exchange plus parameterizations of the short-range 
and intermediate range parts of the potential. The Nijmegen group [15] has 
constructed potentials based on meson exchange and models parameterized 
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in similar ways as the Argonne potentials. Another important difference be- 
tween e.g., the Bonn potentials and the Argonne and Nijmegen potentials is 
the strength of the much debated tensor force [18]. Typically, the Bonn po- 
tentials have a smaller D-state admixture in the deuteron wave function than 
the Argonne and Nijmegen potentials, as well as other potential models. A 
smaller(larger) D-state admixture in the ground state of the deuteron means 
that the tensor force is weaker (stronger). The strength of the tensor force has 
important consequences in calculations of the binding energy for both finite 
nuclei and infinite nuclear matter, see e.g., the discussion in Ref. [19]. A poten- 
tial model with a weak tensor force tends to yield more attraction in a nuclear 
system than a potential with a strong tensor force. The second aim of this 
work, is therefore to compare various nucleon-nucleon potentials, in order to 
study the role played by the tensor force. For this purpose we will work with 
three older versions of the Bonn potentials defined in table A.l of ref. [9]. The 
only essential difference between these potentials is the relative strength of 
tensor force, which make the role of the tensor force more transparent in our 
investigation. They are recognized by the labels A, B and C with the former 
carrying the weakest tensor force. Such a set of potentials, fit to the same set 
of scattering data, allow us therefore to study the role played by the nuclear 
tensor force in a many-body approach like the present. Moreover, these po- 
tentials are constructed in exactly the same way, the only differences being 
the value of various meson variables (coupling constants, masses and cutoff 
energies). One may therefore be able to ascribe possible differences to the role 
played by certain mesons, or certain components of the NN interaction. In 
addition we have also done extensive calculations using the charge-dependent 
version of the Bonn potential models (CD-Bonn), see Ref. [14], and two poten- 
tials of the Nijmegen group [15], one meson-exchange potential (Nijm-OBE) 
and a local potential with a non-local contribution to the central force, i.e. the 
model Nijm-I of Ref. [15]. The reason for these choices is due to the fact that 
these potential models represent the most recent fits to the scattering data. 

The PPHH-diagram nuclear matter calculations are closely related to the 
MBHF calculations, both employing a model space (k < ku)- In section 2 
we will present the MBHF method and the results obtained with it. The ba- 
sic method presented in section 2 is also used in the PPHH-diagram method 
discussed in section 3. The results of the PPHH-diagram calculations for nu- 
clear and neutron matter are presented and discussed in sections 4 and 5, 
respectively. In section 6 we summarize and conclude. 



2 Model Space Brueckner-Hartree-Fock Approach 

The basic formalism of the MBHF has been exposed elsewhere, see e.g., Refs. 
[4,5]. Here we will therefore only briefly sketch the essential ingredients of the 
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method. 



2.1 The Model Space G -matrix and Single-particle Spectrum 



Following the conventional many-body approach, we rewrite the full Hamilto- 
nian H = T + V, T being the kinetic energy and V the bare NN potential, 
in an unperturbed part H Q = T + U and an interacting part Hi = V — U 
such that H — T + V — Hq + Hi, where we have introduced an auxiliary sp 
potential U. If U is chosen such that Hi be small, perturbative many-body 
techniques can presumably be applied. A serious obstacle to any perturbative 
treatment is the fact that the bare NN potential V is very large at small in- 
ternucleon distances, a fact which renders any perturbative approach highly 
prohibitive. To overcome this problem, one introduces the reaction matrix G, 
given by the solution of the integral equation 



v m + £ Vi 



tjmn 



Q M (m,n) 



with the model-space Pauli exclusion operator 



Q M (m,n) 



1, min(/c m , k n ) > kp and 
max(/c m , kn) > kM ■, 

0, otherwise. 



(2) 



The model-space sp spectrum is defined self-consistently by 



_M 



+ I n M 

tk+U k , 



(3) 



where t k is the kinetic energy and the sp potential is 



u 



M 



E h<kF (kh\ G M (u = ef + eff ) \kh) , k < k M , 



M , r .M\ 



(4) 



0, 



k > k M - 



The above MBHF approach is a generalization of the well-known BHF ap- 
proach, the latter corresponding to a special case of the former with ku = kp 
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Fig. 1. The MBHF sp potential U/f shown for different values of ku using the Bonn 
A potential (k F = 1.4 fur 1 ). 



and commonly referred to as "standard" BHF. For the BHF method, only 
the kinetic energy (free particle spectrum) is included in the sp energy for 
momenta greater than hp, while below kp the sp potential u^f is included as 
well. In the MBHF method, the corresponding boundary is at &m, beyond 
which one uses the free-particle spectrum. The BHF and MBHF sp energies 
are constructed so that they are always real. To ensure this property, we need 
to require that the MBHF spectrum e M (k) has a sp potential which is 
always attractive. (Otherwise the propagator in G M will have poles, leading 
to a complex G M .) In Fig. 1, we show the sp potential for some typical 
values of ku- Similar results are obtained using the other potentials. Generally 
speaking, when /cm is small the sp potential is attractive. As ku increases, the 
MBHF sp potential becomes more repulsive. For ku = 3.8 fm - 1 , uf is still 
attractive for small k, but becomes less attractive for larger k and is close to 
zero at k — ku- For fc M > 3.8 fm -1 , the sp potential becomes repulsive. When 
this happens, we will have a crossing between the free particle spectrum and 
the MBHF spectrum. This will in turn give rise to a sp potential which is no 
longer real, unless some ad hoc procedure is imposed such as suppressing its 
imaginary part. To prevent such crossings, we have chosen ku ~ 3.8 fm -1 as 
an upper bound for k M , which defines our model space. This caveat ensures 
that the MBHF G-matrix is always real. 
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The MBHF approach is basically an intermediate scheme, in the sense that 
only two-particle scattering states outside the model space are accounted for. 
This means that both particles must have momenta larger than kp and that 
at least one of the fermions in the two-particle wave function has momentum 
larger than ku- The model space Pauli operator defined in Eq. (2) exclude 
scattering states with both momenta below ku- In the lowest order, MBHF 
calculations for the nuclear matter binding energy is performed in two steps. 
First, Eqs. (1) and (3) are solved self-consistently. In the second step a "stan- 
dard" Brueckner type calculation (see Eq. (5) is performed except that the 
MBHF spectrum defined in Eq. (4) is used for k F < k < k M - 



2.2 MBHF Results 



We have performed lowest order MBHF calculations for various /cm values 
within the upper bound /cm = 3.8 fm -1 . In Fig. 2 we present the average 
energy per particle as a function of ku for three typical Fermi momenta, 
using the Bonn A and C potentials. For small Fermi momenta (kp = 1.2 
fm -1 ) the results for the Bonn A potential are very insensitive to the choice 
of kM- For a wide region of k M (2.0-3.5 fm" 1 ) the energy variation is less than 
0.5 MeV. For larger Fermi momenta we see that the variation is somewhat 
stronger, however, we find a "saturation" behaviour in the vicinity of kM ~ 
3.0/m _1 . This is a fortunate and desirable result. Here we have a minimum 
for E/A and moreover it is here that E/A seems to be least sensitive to k M . 
Thus we believe that /cm = 3.0 fm -1 is an optimum choice for the MBHF 
calculations using the Bonn A potential. Note that intermediate states with 
large momentum components (k > 3 fm -1 are induced by the short-range part 
of the NN interaction [4,21,22]. Thus, using kM ~ 3.0 fm" 1 the intermediate 
states induced by the short-range repulsion are included in the model space G- 
matrix. For the Bonn C potential we get very similar results, although there is 
a slightly larger /c^-dependence than for the A potential. For both potentials 
we find that E/A varies no more than 0.5 MeV for kM between 2.5 fm -1 and 
3.2 fm -1 . Similar results are also found for the Bonn B potential. 

One has to keep in mind that the strengths of the tensor force of the potentials 
discussed are different, as Bonn A contains the weakest tensor force compo- 
nent of the three potentials considered, while Bonn C contains the strongest 
one. Since all three potentials reproduce the same set of scattering data, a 
potential containing a weaker tensor component needs a strong central com- 
ponent. Intermediate states induced by the tensor force are relatively more 
strongly dominated by lower momentum components than states induced by 
the central force [9]. Therefore, one can expect that the energy contribution 
from scattering into states with momentum components lower than kM will 
be larger for potential C than for A. This is also reflected in the larger k M - 
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Fig. 2. Average energy per nucleon as a function of ku for MBHF calculations 
using version A and C of the Bonn potentials. The results are for four typical Fermi 
momenta. 



dependence in the calculations using potential C. The results using the other 
potentials are very similar and we find that it is a good approximation to use 
the model-space boundary ku =3.0 fm _1 for all potentials considered. In Fig. 
3 the saturation points for calculations using ku =3.0 fm _1 are indicated by 
filled symbols while the saturation points from standard BHF calculations, 
which correspond to ku = are indicated by open symbols. We see that 
the MBHF approach yields about 3 MeV additional binding energy compared 
to the BHF approach. The saturation density is slightly increased. The results 
show the well-known fact that within lowest-order Brueckner theory one is not 
able to reproduce the empirical data for nuclear matter. In the next section 
the MBHF method is extended to include additional many-body terms, such 
as the PPHH diagrams. 
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Fig. 3. Nuclear matter saturation points obtained in BHF calculations (opens sym- 
bols). MBHF calculations using ku = 3.0 fm -1 are indicated by filled symbols. The 
results using versions A, B and C of the Bonn potential are shown by squares while 
the results for the charge dependent Bonn potential (CD-Bonn) are indicated by 
circles. The triangles represent the results using Nijmegen potential type I (Nijm-I) 
and the results from the Nijmegen meson exchange version (Nijm-OBE) are denoted 
by triangle+square. 

3 Summation of Particle-particle and Hole-hole Diagrams for Nu- 
clear Matter 



The details of the PPHH-diagram method for nuclear matter can be found 
elsewhere [5]. However, to facilitate the presentation and discussion of our 
present work, we will give a brief review of the method. In perturbative ap- 
proaches to nuclear matter, the well-known BHF theory has been the standard 
starting point. The effective two-body interaction in nuclear matter has then 
been given by the so-called G-matrix, which includes all ladder-type diagrams 
with particle-particle intermediate states to infinite order. However, in terms 
of the G-matrix it is nevertheless only a lowest-order theory. Other many-body 
contributions like screening terms or more complicated many-body terms are 
missing. The ground-state energy shift A£"o in terms of the G-matrix is rep- 
resented by the first-order diagram of Fig. 4 and reads 

AEq HF = Y,n a n b (ab\G BHF (uj = e a + e b )\ab). (5) 

ab 
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In Eq. (5) the n's are the unperturbed Fermi-Dirac distribution functions, 
namely —1 if k < kp and =0 if k > kp where kp is the Fermi momentum. 
The sp energies are denoted by e, and are determined self-consistently using 
the BHF theory. 

Fig. 4. First order contribution to the ground-state energy shift AEq in BHF theory. 



The Brueckner-Hartree-Fock G-matrix G BHF contains repeated interactions 
between a pair of "particle" lines, as illustrated by the diagrams of Fig. 5. Note 
that they are so-called Goldstone diagrams, with an explicit time ordering. The 
third-order diagram (a) of Fig. 5 is given by 

Diag\at) = (-) -. (6) 

2 (e a + e b — e m — e n )(e a + e b — e r — e s ) 



Here m, n, r, s are all particle lines, and V^ki represents the anti-symmetrized 
matrix elements of the NN interaction V. 

To the same order, there is also a diagram with hole-hole interactions, as shown 
by diagram (f3) in Fig. 6, where a, b, r, s are all hole lines. This diagram is not 
included in standard BHF calculations of nuclear matter, for the following 
reason. In earlier times, nuclear-matter calculations were based on, by and 
large, the so-called hole-line-expansion. The essence of the hole-line approach 
is that diagrams with (n+1) hole lines are generally much smaller than those 
with n hole lines. With this criterion, diagram which has 3 (independent) 



o--o * 
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J_ 







(a) 



Fig. 5. Goldstone diagrams contained in the Brueckner-Hartree-Fock G-matrix 

qBHF 



10 




(P) 

Fig. 6. Goldstone diagram with repeated interactions between hole lines. 

hole lines would be negligible compared with diagram (a) which has 2 hole 
lines. Thus the former can be neglected. This approximation has, however, not 
been rigorously checked. To investigate the validity of this criterion, it may be 
useful to actually calculate diagrams like 

A motivation behind the PPHH-diagram method of nuclear matter is also to 
include diagrams with hole-hole correlations like diagram (f3). In so doing, 
the above time-ordered formalism is no longer convenient. We need to use 
the Green's function formalism, where all time variables are integrated over 
the same time interval from — oo to oo. Diagrams (a) and ((3) are both time 
ordered diagrams. In the Green's function formalism, the sum of these two 
diagrams becomes a single diagram, namely diagram (7) of Fig. 7. Note that 
we choose to draw Green's function diagrams differently from the time-ordered 
ones. The diagram of Fig. 7 is given by 




(Y) 

Fig. 7. Green's function diagram discussed in text. The states denoted by r and s 
can either be particles or holes. 
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— 1 00 

Diag.in) = 2ri J d ^° + 



x tr[-F ab (u)V abmn F mn (u)V mnrs F rs (uj)V rsab \, (7) 
where F represents the free PPHH Green's function 



F ( a _ (1 -w a )(l - Wh) n a n b 

ab{ ' co-(e a + e b ) + iO+ co-(e a + e b )-tO+' W 



and V is related to V by 



Vabcd = -V abc d- (9) 



In the above, the symbol tr (trace) means that all indices are summed over. 

The advantage of using this Green's function formalism is that all the PPHH 
diagrams, of the general structure shown by Fig. 7, can be expressed in a 
simple algebraic form. In fact the entire PPHH-diagram series of this type 
becomes 



1 oo 

"PPHH" = ^— J J duje tu)0+ tr[F{u)V 

-oo 

+ l -{F{u)Vf+ l -{F{u)Vf + - 
-1 }d\ 



-oo 

+ {F{oo)V\) 2 + {F{uo)V\f + •••]. (10) 



Here, we have introduced the integration over the parameter A, in order to 
cast the summation over all terms into a simple geometrical series, without 
the factors \, | and so forth. 

The above expression is not yet suitable for nuclear-matter calculations. As 
is well known, the free NN interaction has a very strong repulsive core. The 
above series, expressed in terms of V, is not well defined as each of its terms 
is very large. A standard way to overcome this difficulty is to convert the 
above series into a series in the terms of the G-matrix. Song et al. [5] have 
shown that within the model space approach the PPHH-diagram series can be 
rewritten in terms of G-matrix vertices. For example a third-order G-matrix 
PPHH diagram is given by the diagram of Fig. 7, with each of its dotted-line 
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interactions replaced by a G-matrix interaction. In this way, the contribution 
of the entire PPHH-diagram series to the nuclear-matter ground-state energy 
can be rewritten as 



i 



AE, 



iPPHH 




JdX £ £ Y m (ab, \)Y^(cd, X) 



m<=(A-2) abcdGP 



(ab\G M (cu = E m \cd), 



(11) 



where P denotes the chosen model space, within which all nucleons are re- 
stricted to have momentum less than ku- The G M and e M are the model- 
space G-matrix and sp spectrum defined in Eqs. (1) and (3), respectively. To 
obtain the latter equation, a contour integration closed in the upper tu-plane, 
see Ref. [5] for further details, has been performed. Moreover, since by con- 
struction the model-space BHF G-matrix has no poles for negative energies 
u>, one avoids thereby imaginary contributions to the energy. This is perhaps 
the main advantage of the model-space approach for calculating the PPHH 
diagrams compared with the SCGF method of Dickhoff, Polls and Ramos [1]. 
However due to this approach, as stated previously, there are other observ- 
ables of great interest which one is not able to calculate within the present 
scheme, namely one-body Green's function, spectral functions and momentum 
distributions. 

Note that the above expression is rather similar to the BHF expression of Eq. 
(5). There the weighting factors are n a rib and the G-matrix is evaluated at 
an energy (e a + e b ). Here the weighting factors are YY* and the G-matrix is 
evaluated at energies E m . The transition amplitudes Y and energies E m are 
given by the RPA equation 



For notational economy we will write /i m (ui, A) = /i m . 

It is important to mention that the normalization of the RPA wave functions 
Y is given by 



£K e f + ef)S ijM + (l-m- n j )XG M (u;)}Y m (ef, A) 




(12) 



with the self-consistent condition 



u = fi m (uj, A) = E m . 



(13) 



£(1 -n a - n b )Y m (ab, \)Y* n (ab, A) 



-1 



(14) 



9fJ. m 



ab 



8ul \U=ftm 
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where \i m = /i m (ui,\). For A = 1, E m and Y m of Eq. (12) are both physical 
quantities [20], namely 



E m = {E--Et% (15) 



and 

Y m (aP) = «|a+a+|^ 2 > (16) 



where Eq and ^ represent the energy and wave function of the nuclear 
matter ground state, and similarly E^~ 2 and the m'th state of the 

nuclear matter system with two less nucleons. 

Eq. (12) is a standard RPA secular equation for the PPHH Green's function, 
and as shown its irreducible vertex function is given by the two-body model- 
space G-matrix G M only. This is, however, inadequate in order to describe the 
Green's function, as there are also one-body terms in the PPHH vertex func- 
tion. Such one-body diagrams are important for the PPHH Green's functions 
(see ref. [20] and references quoted therein). To first order in G M , we should 
also include in the vertex function one-body diagrams. The resulting vertex 
function, denoted by T, is then composed of the three diagrams of Fig. 8. Here 
diagram (a) is the two-body G M term, while diagram (b) is the one-body G M 
term. We have added and then subtracted a one-body sp potential u M to the 
Hamiltonian, namely H = T + V = (T + u M ) + (V -u M ) = Hq + Hi. Diagram 
(c) of Fig. 8 is the — u M term. 









+ 




+ 


- - -X 






V 




k ) 


k > 


k i k 



(a) (b) (c) 



Fig. 8. Structure of PPHH vertex function T used to calculate PPHH diagrams in 
this work. The two-body part is the model space G-matrix G M (a). Diagram (b) is 
the one-body term from G M and diagram (c) is the one-body potential introduced 
u M in the Hamiltonian. 

The above concern about the vertex function suggests that for nuclear matter 
calculations we should consider a generalized PPHH-diagram series, as repre- 
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Fig. 9. An example of PPHH diagrams included in our calculations. 



sented by Fig. 9. There each "blob" represents the vertex function T of Fig. 
8. In this way the nuclear-matter PPHH diagrams include not only two-body 
interactions via Gm but also one-body insertions. The all-order sum of these 
generalized PPHH diagrams can still be readily performed, with the final re- 
sult still given by Eq. (9), except that the amplitudes Y and energies E are 
now given by a modified RPA equation, namely 

£[(ef + <;')'W + (1 - th - nj )\r(u;)}Y m (ef, A) 

ef 

= Hm(uj, X)Y m (ij, A). (17) 

A computational subtlety may be mentioned. The above RPA equation are to 
be solved with the self-consistent condition (13). It is rather complicated to 
do so numerically, and it was mainly because of this consideration that earlier 
calculations [5] treated this self consistency by way of a perturbation method. 
In the present work we have treated this self-consistency exactly. 

The inclusion of the one-body insertions for the PPHH diagrams has played 
an important role in our calculation. When such insertions are included, the 
wound-integral effect is present in the PPHH-diagram calculation. Because of 
the interaction among nucleons, particularly the short range correlation among 
them, nucleons in the model space are part of the time outside the model space. 
This "absence" from the model space is represented by the wound integral. 
At low density, the wound integral is generally small and the normalization 
factor in Eq. (14) of the RPA amplitudes is close to -1. As the density increases, 
the wound-integral becomes larger and the normalization can become much 
smaller, such as -0.8. As a result, the PPHH-diagram contribution to the 
nuclear matter potential energy, as given by Eq. (11), is suppressed as density 
increases. In short, the inclusion of the one-body insertions has been essential 
for the determination of the nuclear matter saturation density, as has been 
observed in ref.[5]. 
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As mentioned earlier, the PPHH diagrams are calculated within the model 
space P, where all nucleons are restricted to have momentum up to ku- A 
main concern has been what value for ku one should use. In section 2 we have 
pointed out that k M should not exceed 3.8 fm" 1 . For k M = k F) the PPHH- 
diagram approach reduces to the standard BHF method. Thus the range for 
kM is between kp and 3.8 fm -1 . 



4 Results of PPHH-diagram Calculations 



We have performed PPHH-diagram nuclear matter calculations for a range of 
ku values within the above upper bound. For each ku value, we first calculate 
the MBHF sp spectrum for various Fermi momenta (/c_f). Then we solve the 
RPA equation in Eq. (12) for several values of A. The resulting wave functions 
Y and eigenvalues E are then used to evaluate the PPHH-diagram sum in Eq. 
(11). In Fig. 10 the average energy per particle E/A is shown as a function 
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Fig. 10. The average energy per nucleon as a function of ku for fixed values of k F . 
The results are obtained with the version A of the Bonn potential. 

of ku for four typical Fermi momenta using version A of the Bonn potential. 
The results are quite stable with variations in ku in the vicinity of ku = 
3.8 fm -1 . The corresponding results using the Bonn B and C potentials are 
shown in Figs. 11 and 12, respectively. We see that the results using Bonn A 
and B are very similar, suggesting that ku = 3.8 fm -1 is an optimal choice 
for calculating PPHH diagrams. However, for the Bonn C potential saturation 
occurs at higher kM values. 

To understand this difference it may be instructive to investigate contributions 
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Fig. 11. Results for the version B of the Bonn potential. See Fig. 9 for further 
explanations. 




Fig. 12. Results for the version C of the Bonn potential. See Fig. 9 for further 
explanations. 

from the different partial waves. In Fig. 13 the contributions from the 3 Si+ 3 D! 
channel are shown. Results obtained with the Bonn C potential are shown by 
solid lines. Comparing these curves with those in Fig. 12, we see that most of 
the energy gained when going from a low to a high kM value comes from this 
channel. Note that in this channel the tensor force plays an essential role. 



In the 1 S channel, where the tensor force does not contribute, the energy shift 
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Fig. 13. Contributions to the average energy per particle from different partial waves. 
The results are obtained with the Bonn A and C potentials. See Fig. 10 for further 
explanations. 

is found to be much smaller. In Fig. 13 the 1 So contribution is shown for /c^ = 
1.8 fm -1 For the other Fermi momenta this channel yields similar results. The 
/^M-dependence of the 1 So channel arises mainly in the normalization factors 
of Eq. (14). Thus, the /^/-dependence in the 1 So channel can be explained 
mainly by the tu-dependence of the vertex function T in RPA equation (17). 
Recall that ^-dependent one-body terms are included in T. When the RPA 
equation (17) is decomposed into separate partial wave channels, the one-body 
terms for the various channels are equal. 

For higher order partial waves the (e.g., I > 4) cu-dependence of the G-matrix is 
known to be very weak and the G-matrix differs little from the bare interaction. 
Therefore, the /cM-dependence for all higher order partial waves is mainly due 
to the uj dependent one-body terms and shows up in the normalization factor. 

For comparison we have also included the 3 Si+ 3 Di and 1 S contribution for 
potential A. We see that the 3 Si+ 3 Di channel for this potential has a much 
weaker /c^-dependence. Note that the Bonn A potential has a weaker ten- 
sor force than the Bonn C potential. This indicates that the strength of the 
tensor force is important for the /^-dependence of the calculations. In the 
1 S channel potentials A and C yield similar results. The small difference can 
be explained by the one-body terms of potential C which generally have a 
stronger ^-dependence than those of potential A. This is reflected in different 
normalization factors for the two potentials. 

To complete our results we have included results for the CD-bonn potential 
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Fig. 14. The average energy per nucleon as a function of ku for fixed values of kp. 
Results are obtained with the CD-Bonn, Nijm-OBE and Nijm-I potentials. 

[14] and two versions of the Nijmegen group [15], the meson-exchange potential 
(Nijm-OBE) and a local potential with a non-local contribution to the central 
force, i.e. the model Nijm-I of Ref. [15]. The results are shown in Fig. 14. We 
see that the results are similar to those using the earlier versions of the Bonn 



19 



potential, in the sense that when the strength of the tensor force increases, 
the stability of the results as functions of ku worsens. 



We find that the BHF calculations for Bonn B and Nijm-I yields similar results. 
Note, however, that the tensor force is stronger for the Nijm-I potential. The D- 
state probability of the deuteron reflects the tensor force in the NN-potential. 
The prediction from the Nijm-I potential is 5.66 % while the Bonn B and C 
potentials yields 5.0 % and 5.6 % respectively. Thus, the tensor force for the 
Nijm-I may be more comparable with that of the Bonn C potential. This is also 
reflected in the large PPHH-diagram contribution for the Nijm-I potential. 
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Fig. 15. Nuclear matter saturation curves (PPHH-diagram calculations) and points 
obtained with various potentials. The solid curves are obtained with the CD-Bonn 
potential, while dotted and dashed curves are obtained with the Nijm-OBE and 
Nijm-I respectively. For each potential calculations using ku = 3.0 fm -1 (upper 
curve) and ku = 3.8 fm -1 (lower curve) are presented. For the Bonn A, B and C 
potentials the saturation points (filled star, open star and diamond respectively) 
are shown for the two ku values considered. In addition the MBHF results are 
presented as in Fig. 3. 

In Fig. 15 we show the results of our PPHH-diagram calculation using ku = 3.0 
fm -1 and 3.8 fm -1 . As seen from the figure the saturation density increases 
slightly when ku increases to 3.8 fm -1 . However, compared to the MBHF 
results all our PPHH-diagram calculations yield a lower saturation density. 
Previous PPHH-diagram calculations for nuclear matter were performed using 
ku =3.2 fm -1 [5,13]. We find that our calculations using the same ku value 
yield similar results. Note that in this work the normalization of the RPA 
wave function is performed by the use of Eq. (14). In earlier calculations [5] 
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a perturbation method has been used in order to normalize the RPA wave 
functions. 



5 Calculations for Neutron Matter 



We have found that the results for the PPHH-diagram calculations are very 
sensitive to the strength of the tensor force. In symmetric nuclear matter the 
contribution from the tensor force to the binding energy comes mainly from 
the 3 Si+ 3 D! channel which has isospin equal to zero. For neutron matter, 
however, this channel does not contribute, and therefore the tensor force is not 
so important. To complete our discussion of the PPHH diagrams we therefore 
include results for neutron matter as well. In Fig. 16 we show the average 
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Fig. 16. The average energy per neutron for neutron matter as a function of ku 
for various densities. The MBHF results are shown by dashed lines while the results 
including PPHH diagrams are shown by solid lines. 

energy per nucleon as a function of ku for various densities. In this figure 
we show only results for the Nijm-I potential, since the other potentials gave 
qualitatively similar results. The results from the PPHH-diagram calculation 
are given by solid lines. For the MBHF results, which are shown by dashed 
lines, the energy does not shift more than about 1 MeV if ku increased from its 
lowest value (k M = k F ) to 3.0 fm -1 . We see that the /cju-dependence is very 
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weak in both the MBHF calculations and the calculations including PPHH 
diagrams. Moreover, the contribution from PPHH diagrams are substantially 
reduced compared to the results for symmetric nuclear matter. 

In Fig. 17 the average energy is plotted as a function of density for the three 
potentials. Results including PPHH diagrams are shown. However, we find 
that the contribution from the PPHH diagrams is small. This has important 
consequences for studies of neutron star properties. Many-body terms such as 
hole-hole terms play a smaller role. It is therefore of interest to perform a self- 
consistent calculation of particle-hole diagrams, of importance for long-range 
effects, in order to assess the relation between particle-particle, hole-hole and 
particle-hole terms in neutron matter or /3-stable matter calculations. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
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Fig. 17. The average energy per neutron for neutron matter using various potentials. 
PPHH diagrams are included using ku = 3.8 fm -1 . The results obtained with the 
CD-Bonn potential are shown by the solid line while the results using Nijm-OBE 
and Nijm-I potentials are shown by dotted and dashed lines respectively. 



6 Summary and Conclusion 

A central ingredient in the PPHH-diagram calculation for nuclear matter is its 
model space size ku- How to determine it has been a long standing problem. 
We have carried out extensive nuclear matter calculations, with various choices 
for ku for both the model-space Brueckner-Hartree (MBHF) approach and the 
approach which includes particle-particle and hole-hole (PPHH) diagrams. We 
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have found that for the MBHF approach the nuclear matter binding energy 
and saturation density vary little with }zm for ku in the vicinity 3.0 fm _1 . The 
location of this saturation region for ku is almost independent of the potential 
used. 

For calculations including PPHH diagrams we find a saturation region where 
the nuclear matter binding energy and saturation density vary rather little 
with ku for ku ~ 3.8 fm" 1 . However, the calculations including PPHH di- 
agrams exhibit a stronger £;M-dependence than the MBHF calculations. Our 
results have shown that the strength of the tensor force is important for the 
&M-dependence. Potentials with weak tensor force yield better convergence for 
large ku than those with stronger tensor force. 

The new potentials which have a \ 2 fit P er datum close to 1 generally yield 
similar results as the older versions of the Bonn potentials. For all poten- 
tials considered, the saturation density is substantially reduced when PPHH 
diagrams are included. 

MBHF calculations for neutron matter are found to be insensitive to the choice 
of ku- Moreover, the contribution from PPHH diagrams is found to be small in 
this case. Thus, combining this with our result for symmetric matter, we have 
shown that the contributions from PPHH excitations in a nuclear medium 
depend strongly on the tensor force component of the NN interaction. 

We have carried out our PPHH-diagram calculations using the MBHF sp 
spectrum. It should be of interest to employ other sp spectra in the summation 
of PPHH diagrams for nuclear matter calculations. As an example, it would 
be desirable to employ a sp spectrum which is obtained in a consistent way 
from the PPHH-diagram calculation of the binding energy. As far as we know, 
how to formulate such a sp spectrum is not yet fully understood. Such a 
formulation is presently being studied by us. 
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